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Abstract 

Using precise lattice QCD computations of the baryon spectrum, we present the first direct 
evidence for the presence of contributions to the baryon masses which are non-analytic in the 
light quark masses; contributions which are often denoted chiral logarithms. We isolate the poor 
convergence of SU{3) baryon chiral perturbation theory to the flavor-singlet mass combination. 
The flavor-octet baryon mass splittings, which are corrected by chiral logarithms at next to leading 
order in SU{3) chiral perturbation theory, yield baryon-pion axial coupling constants D,F,C and 
H consistent with QCD values; the first evidence of chiral logarithms in the baryon spectrum. 
The Gell-Mann-Okubo relation, a flavor-27 baryon mass splitting, which is dominated by chiral 
corrections from light quark masses, provides further evidence for the presence of non-analytic light 
quark mass dependence in the baryon spectrum; we simultaneously find the GMO relation to be 
inconsistent with the first few terms in a taylor expansion in nis — mi, which must be valid for 
small values of this SU{3) breaking parameter. Additional, more definitive tests of SU{3) chiral 
perturbation theory will become possible with future, more precise, lattice calculations. 
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I. INTRODUCTION 



Quantum chromodynamics (QCD) is one of the fundamental gauge theories of the stan- 
dard model of particle physics, encoding the interactions amongst quarks and gluons. At 
high energies, the theory exhibits the property of asymptotic freedom where the coupling 
between the quarks and gluons runs to zero as the interactions are probed with larger mo- 
mentum transfer. Conversely, at low energies, at a scale of Aqcd ~ 1 GeV, the coupling 
between the quarks and gluons becomes 0{1), and the theory is no longer amenable to a 
perturbative treatment; the quark and gluon degrees of freedom are bound into the observed 
hadronic degrees of freedom, the protons, neutrons, pions, etc., which leave only subtle clues 
about the underlying fundamental theory of QCD. 

These properties of QCD, as well as many others, are now well established thanks to 
a variety of techniques that have been developed to understand the rich phenomena that 
emerge from the theory. One of the most important tools is lattice QCD, a numerical solution 
to the theory, performed on a discrete, Euclidean space-time lattice. With algorithmic 
advances and ever growing computing power, state of the art lattice QCD calculations are 
performed at several lattice spacings, with moderate physical space-time volumes and with 
dynamical light quark masses at or near their physical values [1-5]. Recently, the ground 
state hadron spectrum, composed of up, down and strange quarks, has been reproduced 
from lattice calculations with a few percent uncertainty [6]. This serves as an important 
benchmark in demonstrating the ability for these numerical calculations to produce precise 
quantitative predictions for hadronic physics observables. Indeed, lattice calculations are 
playing an important role in many areas of both nuclear and high energy physics [7]. 

In addition to this numerical solution to QCD, a variety of analytic methods have been 
developed to understand the low-energy regime of the theory. The most prominent method 
is chiral perturbation theory (xPT) which exploits an approximate global symmetry of 
QCD [8]. For the up, down and strange quarks with masses less than Aqcd, the QCD La- 
grangian is approximately invariant under global chiral transformations of the quark fields 
such that the theory has an approximate SU (3) l x SU (3)/? chiral symmetry, which becomes 
exact in the limit the quarks are massless. This approximate chiral symmetry is sponta- 
neously broken to the SU{3)v subgroup by the QCD vacuum giving rise to the pion octet 
pseudoscalar pseudo-Nambu-Goldstone bosons, the pions, kaons and eta. The realization 
of this chiral symmetry, as well as its spontaneous and explicit symmetry breaking, can 
be described by constructing a chiral Lagrangian which contains this pion octet as well as 
other SU{3)v hadron multiplets as explicit degrees of freedom. In the chiral limit, the pion 
octet become exact Nambu-Goldstone bosons which have only derivative couplings to them- 
selves and other hadrons. This theory, xPT, is non-renormalizable and contains an infinite 
number of operators whose forms are constrained by the global symmetries of QCD, while 
the coefficients of these operators, the low energy constants (LECs), are unconstrained and 
must be determined by comparing with experimental data and/or the results of numerical 
lattice QCD calculations. The quantitative relevance of these operators are dictated by an 
expansion in the soft momentum of the pion octet and the light quark masses suppressed 
by the chiral symmetry breaking scale, A^; hadronic observables can be computed at low 
energies to any fixed precision by keeping operators to a given order in the chiral expansion, 
thus requiring determination of only a finite number of the LECs. 

One of the principle applications of xPT has been to determine the light quark mass 
dependence of various hadronic observables, the simplest of which is the light hadron spec- 
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trum. The motivation comes from the significant numerical cost of performing lattice QCD 
calculations at the physical values of the up and down quark masses. xPT can be used to ex- 
trapolate the numerical lattice QCD results to the physical values of the light quark masses, 
in the process determining some of the LECs associated with the quark mass dependent 
operators. This program has been very successful when applied to the Nambu-Goldstone 
meson spectrum and decay constants, see Ref. [9] for a review, beginning with the first 
significant comparison of lattice QCD results with xPT [10]. 

The comparison with the light baryon spectrum has been wrought with more significant 
difficulties and the overall convergence, and usefulness of the SU (3) baryon xPT is in ques- 
tion. These challenges are not unexpected; first, there is a dense spectrum of low lying 
excited states, introducing new scales in the theory; second, while the expansion parame- 
ters of xPT in the Nambu-Goldstone meson sector are given by ~ ^if,7r,r;/^x' small 
expansion parameter when the baryon fields are included becomes e ~ mK,n,r^/^x i^'^]- 
the physical kaon, e ~ 1/2 and from general expectations of asymptotic series, one does not 
expect this theory to have a controlled perturbative expansion. A few recent comparisons 
of SU (3) baryon xPT to numerical lattice QCD results have led to the conclusion the three 
fiavor chiral expansion is failing to provide a controlled, convergent expansion [12-14] . The 
issues of convergence are not limited to the three fiavor expansion [15] and recent analysis in- 
dicates the range of expansion of the two flavor theory, considering only an expansion about 
the limit of vanishing up and down quark masses, extends only to m^"'^ ~ 300 MeV [12, 16]. 

These challenges have led to a number of efforts to reorganize the expansion for baryon 
xPT. The initial approach is known as heavy baryon xPT (HBxPT) which treats the baryons 
as nearly static fields allowing for an expansion in inverse powers of the baryon mass [17, 18], 
modeled after the heavy quark effective theory [19]. This led to significant phenomenolog- 
ical successes which are partly reviewed in Refs. [20, 21]. Early on, it was recognized the 
convergence of the theory would be problematic because of the large contributions from 
kaon and eta loops in various observables. A new regularization scheme was proposed, the 
introduction of a (chiral symmetry violating) long range regulator, eg. a dipole regulator, to 
soften the contribution from the kaon and eta loops [22]. When applied to chiral extrapola- 
tions of lattice QCD results, this led to some successes in simultaneously describing both the 
numerical results and physical observables [23, 24]. An additional reorganization of the chi- 
ral expansion, equivalent to a resummation of the leading kinetic corrections to the baryon 
propagators was constructed and has become known as infrared regularized baryon xPT [25] ; 
with several offshoots to deal with renormalization of higher loop corrections [26, 27]. Lat- 
tice QCD calculations in the last few years have also made feasible the use of the SU{2) 
expansion for hyperons [28, 29]. In this work, we further examine a new application of an 
old idea: combining the large Nc expansion [30, 31] with the SU{3) chiral expansion [32-36]. 
This approach has a few formal advantages over the other methods. In the large limit, 
there is an extra symmetry, the contracted spin-fiavor symmetry [32, 33], allowing for an 
unambiguous field-theoretic method to include the low lying decuplet baryon resonances in 
the theory; in the large Nc limit, the spin- 1/2 and -3/2 baryons become degenerate and 
infinitely heavy. Further, while the large Nc and SU (3) chiral expansions on their own may 
not provide well converged effective theories, the combined expansions may prove sufficient 
for a controlled perturbative expansion. This approach was first explored in Ref. [37] where 
it was demonstrated the predictions from the combined large Nc and SU (3) expansions on 
the baryon spectrum are well met for a range of light quark masses. 

Having a controlled expansion is necessary but not sufficient to claim success. The princi- 
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pie prediction from xPT are the contributions to hadronic observables which are non-analytic 
in the light quark masses. The masses of Nambu-Goldstone boson is given to leading order 
by the Gell-Mann-Oakes-Renner Relation [38], m^j = B{mi + mj), with a meson composed 
of a quark-anti-quark pair of (anti) flavors i and j and is the mass of a quark with flavor 
i. Therefore, in xPT, the non-analytic light quark mass dependence arises from pion-octet 
loops, which often contribute ln(m|^^^^) terms to hadronic observables, and are commonly 
referred to as chiral logs. These contributions can not arise from a finite number of local 
counterterms but only from the long range contributions from the light pion octet degrees 
of freedom, the pion cloud. Isolating this predicted light quark mass dependence in lattice 
QCD results has been a major challenge for many years. The definitive identification of 
these contributions is hailed as a signal that the up and down (and strange) quarks are 
sufficiently light that the lattice results can be described accurately by xPT. This task has 
proved to be very challenging, as often, these non-analytic light quark mass contributions 
are subleading, or masked by other systematics. 

In this work, we present for the first time, direct evidence of non-analytic light quark mass 
dependence in the baryon spectrum. As will be discussed in this article, this work is not the 
definitive work on the subject, as there are many systematics which must be resolved, but 
this is an important first step in the quest for chiral logs. 



II. THE HEAVY BARYON CHIRAL LAGRANGIAN AND THE LARGE EX- 
PANSION 

A. Heavy Baryon Chiral Lagrangian in the Expansion 

The three flavor heavy baryon chiral Lagrangian at leading order (LO) in the momentum 
expansion and to first order in the chiral-symmetry breaking quark mass matrix M.q = 
diag(m„, md,ms) is given by [17, 18], 

£ = Tr fi, {iv ■V)B,- fl^ {iv ■V)T,^- ^Aq Tr B,B, + ^Aq 7;^T„^ 

+ 2D Tr {B,S^ {A,, B,}) + 2F Tr {B^S^i [A^, B,]) 
+ C {T^A^B, + B,A^T^) +2'HT^S:A,T,^ 
+ 2aB Tr {B,B,) TxM+ - 2gt T^T.^YyM^ 

+ 26^Tr {B, + 26fTr (B„ [M+, B,\) + 2bTT^^M+n^ (1) 

where the spin- 1/2 octet baryon fields B^ and spin-3/2 decuplet baryon fields arc two- 
component velocity-dependent baryon fields which are related to the usual four-component 
relativistic Dirac spin baryon fields B and by 

Tl^{x) = i±ie'^o''-"T^(a;) . (2) 

The mass Mq is the flavor-singlet mass of the baryon octet and decuplet baryons in the 
SU (3) chiral limit nig 0. Specifically, 

Mo=^(M8)-i(Mio), (3) 
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where (Mg) and (Mio) arc the average flavor-singlet masses of the spin-1/2 flavor-octet 
baryons and the spin-3/2 flavor-decuplet baryons, respectively, in the chiral limit. In the 
large expansion, Mq is 0{Nc) for baryons with quarks. The leading heavy baryon 
chiral Lagrangian also contains the flavor-singlet hyperfine mass splitting 

Ao=(Mio)-(M8), (4) 

which is proportional to the total spin-squared of each baryon multiplet. The mass 
parameter Aq is 0{1/Nc) in the l/Nc expansion. The SU{3) flavor representations of the 
QCD baryons are the flavor-octet 

and the completely symmetric rank-3 flavor-decuplet Tijk, normalized such that T^uu = A"^''". 
The heavy baryon chiral Lagrangian also contains four independent baryon-pion couplings, 
the axial couplings D, F, C and "H. The couplings D and F describe the usual baryon- 
octet pion couplings; C describes pion couplings between octet and decuplet baryons; and 
H describes the pion coupling of the decuplet baryons. The pion octet fields 




n = TT^T'^ 



(6) 



appear in the heavy baryon chiral Lagrangian in the nonlinear representation = E = 
g2in//^ where / ~ 130 McV is the pion decay constant in the chiral limit. The vector and 
axial vector pion combinations 

V, = \ + ed,i) , (7) 

appear in the baryon-pion couplings and through the baryon covariant derivative — 
+ iV^. In the heavy baryon chiral Lagrangian, S!^ is the spin operator which acts on the 

spinor portion of the baryon field. 

Additional dependence on the pion field enters through the quark mass matrix spurion 

M+=\ [iM.^i + i^M,i^) . (8) 

In this work, we compare with lattice computations performed with degenerate u and d 
quark masses = ma — mi, so the quark mass matrix reduces to 

= ^ {2mi + m,) U + (m/ - m,) . (9) 

There are two flavor-singlet contributions to the baryon masses with one insertion of the 
quark mass matrix coming from the terms proportional io and gt- There are also three 
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flavor-octet contributions to the baryon masses with a single insertion of the quark mass 
matrix, proportional to 6/), bp and 6^ (called be previously [39]). 

The 1/Nc expansion [30] for baryons [31] leads to the emergence of a spin-flavor sym- 
metry [32, 33, 35] for large- A^c baryons. In Ref. [40], the heavy baryon Lagrangian was 
formulated in the l/N^. expansion. Relations amongst the coefficients in the heavy baryon 
chiral Lagrangian occur at leading and sublcading orders in the 1/Nc expansion, which re- 
duces the number of independent chiral coefficients in the heavy baryon chiral Lagrangian 
at leading and subleading orders in In addition, there exists a planar flavor sym- 

metry [40] at leading order in 1/Nc, which relates flavor-singlet to flavor-octet parameters 
at this order, further reducing the number of independent chiral coefficients in the heavy 
baryon chiral Lagrangian at leading order in the l/N^ expansion. In particular, planar QCD 
flavor symmetry relates the flavor-singlet quark mass parameters ctb and ctt to the flavor- 
octet quark mass parameters bjj, bp and 6^ at leading orders in l/N^.. The flavor-octet and 
flavor-singlet quark mass parameters are given in terms of the coefficients 6„ of the spin-0 
flavor-octet 1/Nc expansion,^ where the subscript n refers to the fact that the corresponding 
operator is an n-body quark operator which is accompanied by an explicit factor of 
jyi-n gj-g^ subleading order in the l/Nc expansion, the mass matrix parameters of the 
heavy baryon chiral Lagrangian for QCD with Nc — 3 are given by 

bD^^b2, bF^^bi + ^b2, bT^-^bi-^b2, 

11 15 
(Tb = + Y^b2 , (Tt = -bi + —b2 . (10) 

The axial couplings D, F, C and H also have an expansion in terms of spin-1 flavor-octet 
coefficients a„ of the 1/Nc expansion. To first subleading order in 1/Nc, the pion-baryon 
couplings of the heavy baryon chiral Lagrangian for QCD with Nc = 3 are related to the 
1/Nc coefficients by [40, 41]^ 

^1 ^11 

-D = -ai , F = -ai + -02 , 

2 6 b 

3 3 

C = — ai , y. — —-ai — -a2. (11) 



B. Mass Relations Ri, R3 and R4 



In Ref. [37], it was argued a better approach to exploring the baryon spectrum was to 
utilize our knowledge of both large Nc as well as SU (3) symmetry which is known to work 
well for the experimental spectrum [43] ; instead of considering the individual baryon masses 
directly, one should explore the light quark mass dependence of various linear combinations 
of the baryon masses, chosen to have definite scaling in terms oi 1/Nc and SU{3) symmetry 
breaking.^ The various linear combinations were determined in Ref. [43]. In Ref. [37], it was 



^ Here, we adopt a simplified notation for the operator coefBcients compared to Ref. [40] . 
^ The 1/Nc operator analysis has recently been extended to the two-body axial current operators [42], such 
as Tr {BA ■ AB) . 

^ Ref. [44] utilized the large relations between operators in baryon xPT to study the baryon spectrum, 
but not the linear combinations constructed to have definite scaling in rus — mi and l/Nc- 
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demonstrated that the predicted scahng with both l/N^ and (m^ — mi) was clearly visible 
in the lattice data. The first few mass combinations had statistically meaningful values over 
the range of quark masses, but there were not enough statistics to resolve all of them. In 
this work, we focus our attention on three of these mass relations, Ri^Rs and i?4.^ These 
mass relations are given by 

Ri = (12) 
Z^j Fiji 

where 

Ml = ^ cijMj = 25(2Mjv + Ma + 3Ms + 2Ms) - 4(4Ma + 3Ms* + 2Ms. + Mo) , 
j 

M3 = ^ csjMj = 5(6Mjv + Ma - 3Ms - 4Mh) - 2(2Ma - Ms* - Mq) , 
j 

M4 = ^ C4jMj = Mat + Ma - 3Ms + Ms , (13) 
j 

and for example R4 — M4/6. 

These relations are designed to isolate various operators in the combined l/N^ and SU{3)- 
breaking expansions. At 0{mq), only relations R1-R4 are non- vanishing. For this reason, 
the relations R5 - Rs are particularly interesting to use with light quark mass extrapolations, 
as the leading contribution begins with the chiral loops at 0{ml^'^). However, even more 
precise results of the baryon spectrum than exist are needed for these relations. Using the 
large Nc expansions through second non-trivial order, and working through next-to-leading 
order (NLO) in the chiral expansion, the relation Ri is given by 

- Ri{mi, rus) = Mo - l-bi + —62 j {2mi + rus) 

_ ^3g^2 _ 5^2^ ^ 3J^{m^,0, iJ,) + 4J^{mK, 0, /x) + J'(m^, 0, //) 



1 2 
-12"^ 



50 



8(47r/)2 

3J^(m^, A, jj.) + AJ'{mK, A, ^) + J^(m^, A, jj.) 



8(47r/)2 

3J^(m7r, -A, jj,) + AJ'{mK, -A, /x) + J^(m^, -A, /j,) 



8(47r/)2 

The non- analytic function is defined as 

J-(m, A, //) = (A^ - + lef In . "^..^ „ = 

V A — V A^ — m^ + iej 



(14) 



-2^- ^n^)"^ ^H^)- ^''^ 



The relation R2 gives at leading order the hyper-fine splitting Aq. For the current lattice data set, this 
quantity provides no further information over the use of Ri . 
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which has the hmits and properties 



J'(0,A,/x) = 

J-'(m, 0, /i) = Tim^ 

-^■(m, A, /i) + 2i7r(A2 - m^)^/^ m < | A| 
-J'(m,A,/x) + 27r(m2- A2)3/2^ m > |A| 



(16) 



For the baryon spectrum, the leading non-analytic light quark mass dependence is encoded 
in this function. As such, it is of particular interest to find evidence of this behavior in the 
spectrum. 

The mass relations R3 and R4 vanish in both the SU (3) chiral and vector limits, making 
them more sensitive to the NLO non-analytic light quark mass dependence. At NLO in the 
chiral expansion, and to the first two non-trivial orders in the large A'^c expansion, these 
relations are given by 

39 Rsimi, ms) = 20 hi (m^ - mi) 

_ 1^20a^ - 5a^) ( Q' ~ "^J^i^K, 0, /^) - ^(m^, 0, yu) 



35 



V (47r/)2 
3J^{m„, A, 11) - 2J^{mK, A, /x) - 7"(m^, A, /i) 
(47r/)2 

3J"(m^, -A, /x) - 2J^{mK, -A, /x) - J"(m^, -A, /x) 

(47r/)2 



[17) 



— 62 {ms - mi) 

,1/2,. , 2\ /3J'(m^,0,//) -2J'(mir,0,/i) - J'(m^,0,/x)\ 
+ Aaia, + a,) ^ — j 

2, 2 f 3J^(m^, A, /i) — 2J^{mK, A, //) — J-'(m^, A, /x) 

In addition to these three mass relations, we also explore the Gell-Mann-Okubo relation 

Agmo = + ^Me - ^Mn - ^Mh . (19) 

Since the quark mass operator contains pieces which transform as both an 8 as well as a 1 
under SU{3) transformations, Eq. (9), there are non- vanishing contributions to the GMO 
relation. However, mass operators which transform as an 8 make vanishing contributions 
to Eq. (19). The leading mass operator which makes a non-zero contribution to the GMO 
relation transforms as a fiavor-27. These corrections can arise either from chiral loops 
or from a mass operator containing two or more quark mass insertions. This makes the 
GMO relation particularly interesting to explore with lattice QCD calculations; the leading 
contribution to this mass relation comes from chiral loop effects which are non-analytic in 
the light quark masses. Experimentally, the GMO relation is found to be 

ASo = 6.45MeV. (20) 



R^imums) = - 
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Each baryon mass in the relation receives non-analytic mass corrections which scale as 
SMb oc N(.m'i^'^ . These large corrections may lead to the expectation that the GMO rela- 
tion receives large contributions from the loop corrections. However, one can show these 
NcTr^J"^ terms are proportional to 1 under SU (3) transformations. Additionally, the rr^J"^ 
contributions transform as an 8 while the ml^^/Nc corrections transform as a flavor-27. This 
provides an extra on top of the chiral suppression, explaining the relatively small value 
of the GMO relation [40]. 

At ncxt-to- leading order in the chiral and large Nc expansions, the Gell-Mann-Okubo 
relation is 



A NLO 
^GMO 



36(47r/)2 



+ 



4aia2 + Oo 



36(47r/) 



J'(m7r, 0, //) - 4:J^{mK, 0, //) + 3J'(m^, 0, //) 

+ 2J'(m7r, A, jj.) - 8J'{mK, A, /j.) + 6J^(m^, A, jj.) 



(21) 



In this article, we will also be interested in the next-to-next-to-leading order (NNLO) 
formula. This can be determined from Ref. [45, 46]. Retaining the subleading Nc relations 
for the quark mass operators, Eq. (10), but only the leading relations for the axial couplings 
(a2 0), the NNLO contributions to the GMO formula are 

aNnlo_ 1 (6r + &f)K-^0' , 23aj fmt-4m% + 3mX 



^GMO 



48 

+ 
+ 



47r/ 



+ 



(- 



V967r/ 192 Mo 



384 Mo 

4 1„/'^2 /,,2 



hi 



13a? \ [ml ln(m^//x2) - 4m|^ \n{m\/ir) + 3m^ hiimy i^) 
[j{m^, A, /x) - AJ{mK, A, /x) + 3 J(m^, A, /x) 



2m; -|- m, 



12 



+ 62 



—9ms + 16m; 
24 



[j{m„, A, /x) - AJ{mK, A, /x) + 3J^(m^, A, /x) 



^(ms - mi) (^Jirrir,, A, //) - J{mK, A, //)) 



- ^(rus - mi) ( m^ln 



m 



K 



— m^ In 



m: 



+ 



16m; -|- 5m<j 
72 



m^ — Am^ + 3m„ ) — 



rris - mi 



m^-m^ 



(22) 



where the function J{m, A, /^) encodes additional non-analytic dependence on the light 
quark masses 



J(m, A,//) = 2AVA2 -m2 + ieln 



A + \/A2 -m^ + ie 
A - VA2 -m^^ie 

+ m^ In 



m 



2A2ln 



4A2 



m^ 



• (23) 
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and has the hmits and properties 



J^(0,A,/i) = 
J{m, 0, /x) = In 



m 



Jim - A = / ^' ^) + 4i7rAVA2 - m^, m < | A| 

^ ' '^^ \ J{m, A, //) - 47rAVm2 - A^, m > A 

III. EVIDENCE FOR NON-ANALYTIC LIGHT QUARK MASS DEPENDENCE 
IN BARYON SPECTRUM 

A. Details of the lattice results 

For this work, the numerical results of Ref. [12] are utilized, which are not the most 
recent but are still the most statistically precise data set available. The lattice calculation 
was performed with a mixed-action composed of domain-wall fermion [47-51] propagators 
generated on the Uf — 2 + 1 asqtad-improved [52, 53], rooted, staggered sea quark configu- 
rations generated by the MILC Collaboration [54] . This particular mixed-action set up has 
been used quite extensively by the LHP [12, 55-59] and NPLQCD [14, 60-74] CoUaboraitons 
as well as some independent works [75-79]. The mixed- action effective field theory, which 
encodes the discretization effects specific to this particular mixed-action, has also been thor- 
oughly developed [80-93]. However, the baryon spectrum results we use in this work exist 
at only a single lattice spacing. There is also reason to believe the discretization systemat- 
ics are small [12, 88, 91] and to the order we are working in the mixed-action EFT, they 
are subleading. For these reasons, we only use the continuum xPT extrapolation formula, 
presented in the previous section. 

To set the scale, we use the latest scale setting by the MILC Collaboration [2] , as detailed 
in Ref. [74]; we first convert the numerical results of Ref. [12] into ri units^ and then use the 
MILC determination of ri(m^^^, m^^^) to convert to physical units. Finally, we perform the 
extrapolations as functions of the quark masses. The quark masses are not renormalization 
scheme or scale independent. However, at a fixed lattice spacing, we can absorb the quark 
mass renormalization into the quantity B, where at leading order, the Nambu-Goldstone 
boson masses are given by mf^ = B{'mi + rrij). We then define lattice quark masses, in 
physical units by 

rf ^mf * = ^(am, + am''') , (25) 

where m*"^^ is the residual chiral symmetry breaking present with the Domain- Wall lattice 
action at finite fifth dimensional extent [94]. We collect all these numerical values in Table I. 

To extrapolate the lattice results to the physical point, NLO xPT [95] is used to determine 
the values of m^"** which reproduce 

^ = 138 MeV , m^^^ = 496 MeV . (26) 



The length scale ri is determined with the heavy quark potential, defined such that rfFin) = —1. 
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TABLE I: The numerical lattice data from Ref. [12], converted to physical units. The uncertainties 
are the statistical and systematic uncertainties from Ref. [12] combined in quadrature. The strange 
quark mass is fixed to the largest am,„ value in all calculations: amf^°- = 0.050; am^'^^ = 0.081. 



am' 
am 
am 



/3 



I 

val 
I 

res 



6.76 
0.007 
0.0081 
0.00160(3) 



6.76 
0.010 
0.0138 
0.00157(1) 



6.79 
0.020 
0.0313 
0.00123(1) 



6.81 
0.030 

0.0478 
0.00101(1) 



6.83 
0.040 
0.0644 
0.00083(2) 



6.85 
0.050 
0.081 
0.00073(3) 



m 



latt 



[MeV] 
[MeV] 
mx [MeV] 



q 



16.8 
320(2) 
640(2) 



26.6 

389(2) 
659(2) 



58.0 
557(1) 
726(1) 



685(2) 
787(2) 



121 

805(4) 
852(4) 



155 
905(2) 
905(2) 



fi?i [MeV] 
i?3 [MeV] 
i?4 [MeV] 
Agmo [MeV] 



1285(6) 
-113(3) 
-39(2) 
5.6(2.3) 



1315(6) 
-100(2) 
-33(1) 
1.8(1.2) 



1454(6) 
-64(1) 
-19(1) 
0.18(48) 



1556(12) 
-41(1) 
-11(1) 

0.13(35) 



1698(13) 

-19(1) 
-4.4(0.6) 
0.13(0.09) 



1769(9) 







FIG. 1: The results of SU{3) xPT analysis of the pion and kaon masses. The left plot displays 
the fit result with the data and the right plot shows the size of the NLO contributions compared to 
LO. This particular fit is performed to just the lightest three values of the light quark mass, and it 
is clear from the plot that the fits including all six points yields consistent restuls. 



It is interesting to note that despite ignoring the issues of quark mass renormalization, this 
yields the values 

mi J*, = 3.0(2) MeV , m}:^ = 99(5) MeV , (27) 

which are remarkably similar to the proper lattice determination of the light and strange 
quark masses [9]. The NLO xPT formula provide a controlled and convergent description 
of both m^r and rriK over the full range of quark masses used, see Fig. 1. 
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B. Large Nc and Consistency of Hyperon Axial Charges 



One of the major failings in the apphcation of SU (3) heavy baryon xPT is a lack of con- 
sistency between the determination of the axial coupling constants, D, F, C and H when 
determined from the baryon spectrum [12, 13] versus a direct calculation of the hyperon 
axial charges [75]. The direct lattice determination^ yields values consistent with the phe- 
nomenological values [96] , while the indirect determination from the baryon spectrum yields 
values consistent with zero. The small values of the axial couplings returned indicate the 
numerical results do not support evidence for the leading non-analytic light quark mass de- 
pendence predicted in the spectrum. This problem is not unique to the S'f/(3) heavy baryon 
XPT extrapolations, with large contributions from kaon and eta loops, but also observed in 
the SU{2) extrapolation of just the nucleon mass. As demonstrated in Refs. [12, 16], for 
m^r > 300 MeV, there are large cancellations between the LO, NLO and NNLO contribu- 
tions to the nucleon mass; in order to accommodate the large negative mass contribution 
occurring at NLO, the leading non-analytic light quark mass dependence, there must be a 
compensating large but positive contribution from the LO and NNLO terms, signaling a 
breakdown of the perturbative expansion for these heavier pion masses. A similar and more 
severe situation occurs for the SU (3) chiral expansion. As we shall demonstrate in the next 
sections, it may be we are simply asking the wrong questions of heavy baryon xPT. 

1. Mass relation Ri 

The mass relation Ri is a flavor singlet mass combination designed to isolate the Mq 
contribution plus higher order chiral corrections. Starting with the lightest quark mass and 
including successively heavier values of rUq, there are four possible ranges of light quark 
masses which can be used to perform the chiral extrapolation analysis. Both LO and NLO 
analyses are performed over all these ranges of light quark masses. In the NLO analysis, the 
subleading in N^. axial coefficient is set to zero, 02 = 0. Several choices of the parameter / 
are taken to explore systematics from higher orders in the chiral expansion: 





(28) 



Prom the LO analysis, the following LECs are obtained 



Mo[LO] = 903(20) MeV, 



bi + ^h] [LO] 



2.73(11). 



(29) 



Extrapolating to the physical values of the light and strange quark masses gives 




(30) 



^ In Ref. [75], gA, <?ee and g~,s were computed which were used to infer the values of D and F. 
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FIG. 2: Representative fits to Ri from LO (left) and NLO (right) HBxPT analysis. The blue 
star is the physical value, not used in the analysis. The upper error band results from a fit to the 
lightest four numerical data and the lower bad is the result extrapolated to the physical value of 



m 



latt 



. Eq. (27). 



which is to be compared with ^R^^^ 
are determined to be 



1093 MeV. Performing the NLO analysis, the LECs 



Mo [NLO] = 899(40) MeV^ 



bi H 62 
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[NLO] = -3.26(70) , ai[NLO] = 0.24(30) 



(31) 



with a determination 



-i?i[NLO] = 1127(50) MeV. 



(32) 



In Fig. 2, representative fits of Ri from LO and NLO are displayed. One may take comfort 
in the consistent values of the LECs Mq and bi + ^^62 between the LO and NLO analyses. 
However, this is not surprising given the small value of oi determined in the NLO analysis. 
This small value is consistent with no contributions from the NLO terms and inconsistent 
with the known phenomenological determination of the axial coupling. This is not surprising 
given the convergence issues observed in the SU (2) extrapolation of the nucleon mass [12, 16]. 
One is left to conclude that the SU{3) heavy baryon xPT does not provide a controlled, 
convergent expansion for the mass combination Ri for the range of quark masses used in 
this work. 



2. Mass relations R3 and R4 

The relations R3 and -R4 both receive leading contributions from flavor-octet mass op- 
erators, vanishing in both the SU{3) vector as well as SU{3) chiral limits. From these 
symmetries, the relations R^ and R4 are more sensitive to the non-analytic light quark mass 
dependence occurring at NLO in the chiral expansion. As with the analysis of Ri, three 
choices of the parameter / are taken to estimate higher order effects, Eq. (28). The LO 
expressions for R3 and i?4, Eqs. (17) and (18) with = 0, do not describe the numerical 
results well; it is clear higher order contributions are necessary for extrapolations of this 
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FIG. 3: The LO and NLO contributions to (left) and R4 (right). A (blue) star is used to denote 
the physical values, not included in the analysis. The particular fit displayed is a combined analysis 
of i?3 and i?4 to the data at the lightest three values of m|"** . 



data. At NLO, the analysis of and R4 becomes correlated. The full covariance matrix is 
constructed as described in Ref. [37]. The numerical results of Ref. [12] are insufficient to 
constrain both the leading and subleading axial coefficients, and so the analysis is restricted 
to the set of LECs 

A = (61,62,01) , (33) 
with a2 = 0. From the NLO analysis, the LECs are determined to be 

61 [NLO] = -6.6(5) , 62 [NLO] = 4.3(4) , ai[NLO] = 1.4(1) . (34) 

Using the leading large Nc relations with 02 = in Eq. (11), this corresponds to 

D = 0.70(5), F = 0.47(3), C = -1.4(1) , H = -2.1(2) . (35) 

The significance of this is prominent; the large value of the axial coupling is strong evidence 
for the presence of the non-analytic light quark mass dependence in these mass relations. 
Further, this is the first time an analysis of the baryon spectrum has returned values of the 
axial couplings consistent with phenomenology.^ 

However, caution is in order. Examining the resulting contributions to and i?4 from LO 
and NLO separately, one observes a delicate cancellation between the different contributions, 
see Fig. 3. Further studies are needed with more numerical data sufficient to also constrain 
the parameter 02 as well as the NNLO contributions. 



C. Gell-Mann Okubo Relation 

The leading contribution to the Gell-Mann-Okubo relation is from a flavor-27, which in 
HB^PT come from the leading non-analytic light quark mass dependence, Eq. (21). For this 



^ Finding values of the axial couplings consistent with phenomenology has not just been a challenge for 
lattice QCD, but also observed in large Nc xPT analysis of the experimentally measured baryon magnetic 
moments [97, 98]. 
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FIG. 4: GMO mass splitting plotted as a function o/mj^**. The * is the PDG point, not included 
in the analysis. The various fits are described in the text. In a given plot, the filled (blue) circles 
denote results included in the analysis while the open (gray) boxes are excluded. 



reason, it is a particularly interesting mass relation to study, as has been done if Refs. [12, 62] . 
In this article, the analysis is taken a step further. Close to the SU (3) vector limit, the GMO 
relation can be described by a taylor expansion in — m;, 

^GMO = ^2 {riis - mif + 4 {nis - mif H (36) 

The leading term proportional to (m^ — mi) must vanish as it transforms as a fiavor-8. 
The {ttIs — niiY contribution is equivalent to an NNLO contribution from HB^PT and the 
{nis — niiY contribution is equivalent to an NNNNLO HB^PT contribution. We demon- 
strate the first few non-vanishing terms in this taylor expansion are inconsistent with the 
numerical lattice data. We further demonstrate the NNLO HB^PT formula can naturally 
accommodate the strong light quark mass dependence, which is dominated by the non- 
analytic contributions. In Fig. 4, four plots are displayed. The first plot (upper left) is the 
result of an NLO analysis of the GMO formula, allowing the axial coupling to be determined 
from the data, resulting in a small, but non-zero value for Oi. The second plot (upper right) 
displays the predicted value of the GMO relation from NLO taking the determination of ai 
from the analysis of -R3 and R^. The third plot (bottom left) shows the result of a taylor 
expansion about the SU (3) vector limit fitting the first two non- vanishing terms. Finally, 
the NNLO analysis is displayed, using Eqs. (21) and (22) with 02 = 0, and taking ai from 
the determination from i?3 and i?4 (bottom right). Only the NNLO analysis is consistent 
with the values of the numerical data over the full range of light quark masses, in particular, 
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the steep rise observed as rrii — > 0, as well as the value of the axial coupling ai determined 
from phenomenology. This is further evidence for the presence of non-analytic light quark 
mass dependence in the baryon spectrum. 

IV. CONCLUSIONS 

In this article, we have presented the first substantial evidence for non-analytic light 

quark mass contributions to the baryon spectrum. This was achieved by comparing the 
predictions from heavy baryon xPT, combined with the large expansion to relatively high 
statistics lattice computations of the octet and decuplet spectrum. The numerical results 
available [12] allowed for a detailed comparison of the mass relations Ri, R3 and i?4 [43] as 
well as the Gell-Mann-Okubo relation. It was demonstrated the poor convergence of SU{3) 
heavy baryon xPT is isolated in the mass relation Ri. An analysis of mass relations R3 and 
i?4 provided for the first time, values of the axial couplings which are consistent with the 
phenomenological determination, signaling significant contributions from non- analytic light 
quark mass dependence in R3 and R4. At leading order in the large Nc expansion, it was 
found 

D = 0.70(5), F = 0.47(3), C = -1.4(1) , = -2.1(2) . 

It was further demonstrated that the Gell-Mann-Okubo relation is inconsistent with the first 
two non- vanishing terms in a taylor expansion about the SU{3) vector limit, and that the 
steep rise in the numerical data, observed as m["** 0, can only be described by the NNLO 
heavy baryon xPT formula which is dominated by chiral loop contributions. Taken together, 
these observations indicate the first significant evidence for the presence of non-analytic light 
quark mass dependence in the baryon spectrum. 

This is not the definitive work however. There arc several known systematics which were 
not addressed in the present article, and require future, more precise lattice results: 

• the numerical data used [12] exist at only a single lattice spacing, 

• a continuum xPT analysis was performed, 

• there may be contamination from finite volume effects [99], 

• the convergence issues need further examination, 

• more precise numerical results are needed to explore mass relations R^ - Rg which 
should be more sensitive to non-analytic light quark mass dependence, 

• results with smaller values of the light quark mass are desireable, 

• the strange quark mass used in this work is known to be ~25% to large [100]. 

Addressing these systematics is beyond the scope of this work. However, current lattice 
calculations underway should be able to explore these mass relations in more detail. In 
particular, the new strategy presented in Ref. [101, 102], where the sum of the quark masses 
is held fixed, ruu + md + ms = fh, for a range of light and strange quark masses, proves very 
promising for comparing with predictions from xPT. Further, the strategy is not limited 
to the spectrum, with similar relations having been recently determined for the baryon 
magnetic moments [103]. 
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